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Abstract 

We provide an algorithm to generate vertices for the Schrodinger functional with an 
abelian background gauge field. The background field has a non-trivial color structure, 
therefore we mainly focus on a manipulation of the color matrix part. We propose how to 
implement the algorithm especially in python code. By using python outputs produced by 
the code, we also show how to write a numerical expression of vertices in the time-momentum 
as well as the coordinate space into a Feynman diagram calculation code. As examples of the 
applications of the algorithm, we provide some one- loop results, ratios of the A parameters 
between the plaquette gauge action and the improved gauge actions composed from six- 
link loops (rectangular, chair and parallelogram), the determination of the 0(a) boundary 
counter term to this order, and the perturbative cutoff effects of the step scaling function of 
the Schrodinger functional coupling constant. 

1 Introduction 

It is well known that vertices in lattice gauge theory are quite complicated especially for the 
pure gauge sector. This is because, on the lattice, while gauge invariance is preserved exactly at 
a finite lattice spacing, the lattice regularization itself breaks the Lorentz symmetry explicitly. 
This is the main difficulty of the lattice perturbation theory. To reduce the risk of errors and to 
alleviate the tedious task of deriving the vertices, it is desirable to have an automatic method. 
A first attempt was made by Liischer and Weisz about twenty yeas ago [Ij. They worked in 
momentum space and performed some calculations by using their algorithm which is restricted to 
closed gauge loops but sufficient for pure gauge theory. About two decades later a new algorithm, 
which we call bottom up algorithm, was proposed by Hart et. al. [2j. A crucial point in this 
generalization is that it can deal with any parallel transporter, not only with closed loops. This 
algorithm also allows us to extend to fermion actions or even smeared HQET actions. In fact, 
it is applied to the NRQCD [3j. A relevant assumption of this algorithm has been translation 
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invariance. Needless to say their algorithm also assumes that the background field is set to a 
trivial one. 

Our main concern in this paper is to extend the bottom up algorithm to the Schrodinger 
functional (SF) [4J, where translation invariance for the time direction is broken, and there is 
an abelian color background gauge field. An essential question here is how to deal with the 
non-trivial color matrix which involves the complicated background field. Actually this issue 
can be solved by noting two important properties of the background field, that is, abelian and a 
linear time dependence. Although the latter property can be relaxed, exploiting this is greatly 
helpful to gain speed up of the calculations. Our main purpose in this paper is an exposition 
of the solution. The procedure to generate vertices consists of two parts, first, generating lists 
which contain information of the non-zero elements of the vertices, and second, writing down 
a numerical expression of the vertex by making use of the lists as an input. In this paper, we 
will describe how to write down the vertex in the Feynman diagram calculation code not only 
for the time-momentum space but also the position space. In the two-loop calculation of the 
SF coupling [5], it was emphasized that the position space calculation has the advantage of 
substantially reducing the computational efforts over the time-momentum space. This is the 
reason why we consider the position space vertex here. As it is the case for the original bottom 
up algorithm, our generalization to the SF is in principle applicable to both fermion and gauge 
actions. 

The rest of this paper is organized as follows. First we review the bottom up algorithm in 
the next section. In section [3l after a brief reminder of the SF, we provide a nice representation 
of link variables and key equations derived from a special property of the background field, and 
then we show the algorithm to generate lists for the SF. This section is the main part of the 
paper. Then, an explanation of how to write down the vertices to a program which calculates 
Feynman diagrams (in particular the Big Mac diagram [5] ) is given in section HI As mentioned 
before we provide it in both the time-momentum and the position space. As an application, in 
section [5l we give a one- loop calculation of the SF coupling for various gauge actions including 
six-link loops. Our conclusion is given in the last section, together with comments on possible 
future work. Some technical details are given in appendices, like how to partially symmetrize 
the vertex for the gauge action in section |El further reduction of the number of lists in |Fl and 
how to obtain the r] derivative of the vertex in \G\ respectively. 

2 Bottom up Algorithm 

In this section, we summarize the position space version of the bottom up algorithm on the 
usual translation invariant lattice (non SF). The explanation here is mainly based on [3j apart 
from the position space. We set up the theory on a hyper-cubic euclidean lattice A with spacing 
a and size x T. We assume periodic boundary conditions for all directions. A coordinate 
X £ A has components with // = 1,2,3,4. In the following we denote the link variable (still 
without background field) by 

U{x,n) = exp{agoq^{x)), (1) 
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and the anti-hermitian gauge fluctuation field by q^{x) = X]a9M(-^)^°' where T"" are anti- 
hermitian and are elements of the Lie algebra of SU{3). They are explicitly given in Appendix 

\M 

A first important point for the automatic operation is how to represent the vertices in a 
program. As an example, we consider a parallel transporter P{C, Xs, x^) along a curve C on the 
lattice. We denote with Xg and Xe the start and end point of the parallel transporter. Taylor 
expansion of the parallel transporter is written as 

oo / Nj, r 

p(/:,..,xe) = 1 + E E K-zU^^ (2) 

r=l ai,--- ,ar ai,--- ,ar j=l 

where a = (//, x) is a combined index labelling links. The (unsymmetrized) vertex V^^'.'.'^^, has 
been factorizeqj into the color factor 

and the reduced vertex, Yaj^...a^. The color factor is independent of the shape of the parallel 
transporter, therefore we consider only the reduced vertex. 

For a usual parallel transporter, almost all elements of the corresponding reduced vertex are 
zero. Therefore usually only non-zero elements are stored in a memory of a machine. This is eas- 
ily realized as follows. First one labels the non-zero elements of the reduced vertex with number 
i properly {i = l,...,Nr), and then one stores its indices, ai[z],--- ,arHj and a corresponding 
non-zero value of the Y, which is represented as the amplitude f[i] in a list 

L('")[i] = {(ai[i],--- ,ar[i\),Xs[i\,Xc[i]; f[i]}. (5) 

The number of lists Nr depends on the expansion order r, and of course on the shape of the 
parallel transporter which one considers. At this point the reader might think Xg and Xe are 
redundant but as we will see these are important when dealing with the multiplication of lists. 
A set of lists of order r, 5^^^ = {-L(^)[i]|i = 1, A^^}, is "equivalent to" the r'th order reduced 
vertex, 

Nr 

E • • • C = E C w • • • Ch/W' (6) 

ai,...,ar i=l 

for a fixed color configuration {oi, • • • , Or-}. Then the parallel transporter in eq.([2|) is re-written 
as 

Pic,x,,x^) = i+J2-^ E c--ECw-Ch/w- (7) 

r=l a\,--- ,ar i=l 

This expression is useful when deriving the multiplication algorithm as we shall see soon. A selH 
S = {5^1 r = 0,1, ...,rinax}) which represents the reduced vertex up to the order r^nax for the 
parallel transporter, is now encoded in a program. 



^ For the SF, however, this nice structure is rendered more comphcated by the presence of the background 
field as we will see in the next section. 
^S'^'') is a subset of S. 
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As an explicit example for the list, let us consider the case of a single link variable. From 
the expanded form 

u{x,^,) = l+f2^-^ c'^--^c---e, (8) 

r=l ai,---,ar 

a corresponding list reads 

L^^^ = {(a,- ■ ■ ,a),x,x + afL]l}, a = (^,x). (9) 

^ V ' 

r terms 

It consists of only one line {Nr = 1) and is the elementary building block of the algorithm, that 
is, something like an initial condition. 

So far we have defined the fundamental elements (the list and set) on which the algorithm 
operates. Next we discuss how to generate the set S for any parallel transporter. First we 
assume that two sets S{C^"'\xi^\x'i^^) (n = 1,2) for parallel transporters xi"^) are 

known, and their elements are given by 

LM[i](£("),4"),xW) = {(aS")H,... ,a(")[f]),xW,xi");/W[^]}. (10) 

Then, we consider a product of the parallel transporters 

P(£« * £(2) ^41)^42) ^ P(£«,x«,x«)P(£(2),42) ^^^^(2) (11) 

where A = shift vector and is chosen such that the end point of the n = 1 parallel 

transporter and the starting point of the n = 2 one are identical. The shift for the n = 2 parallel 
transporter in eq. pij) is taken into account for the list structure as 

L^^M{6'\xf\xfy) lMW(£(2),x(2)+A,42)+a) 

= {(af)[i],--- ,df[i]),42)+A,x(2)+A;/(2)[i]}, (12) 

where all x are shifted by A and a = {fi, x + A). The question now is how to obtain the set for 
the product of the parallel transporter in eq. (|ll|) . 

5(£«,x«,4^))*5(£(^),xF,x(2))^S(£«*/:(2),x«,4^)+A), (13) 

from the two given sets S {C^"'\ xi"'\ x^^ ) with n = 1,2. Since any parallel transporter is 
composed of the elementary single link variables, by starting from a set for the link variable 
whose elements are shown in eq.Q and by repeating the above multiplication one can obtain a 
set for arbitrary parallel transporters. This is the origin of the name "bottom up" . The algorithm 
can be understood by looking at the actual multiplication of two parallel transporters. The right 
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hand side of eq.lfTT]) is expanded as 



1=1 1=1 



X 

s 



r=l ai,---,ar k=l 

After the first equal-sign we inserted eq.([7]). In the last step, we have used the fact that the 
color factor is independent of the shape of the parallel transporter and have combined the three 
summations (over s, i and j) into that over k. In short, we did a relabelling of the indices, 
and rewrote the amplitude factor. Finally the resulting lists L^^^[k]{£^^^ *C^'^\ are 
created by putting the new label structure and the new amplitude. The algorithm is summarized 
as follows 

• Relabelling: 

{a'i^i],--- ,a''P\i],af\j],--- ,a^i^Ji]} — >{ai[k],--- ,as[k],as+i[k],- ■ ■ ,ar[k]}, (15) 

• Amplitude part: 



r! 



s!(r — s)\ 

Creating list: 



f^'^\i]f^'^[3]^m, (16) 



l(^)[^](£«,x«,x«) *l(^-)[j](£(2),42) ^^^^(2) _^ 

lM[A:](£« * £(2)^ ^(1)^^(2) ^ ^) ^ . . . ,«^[A;]),4i),42) + A; f[k]}. (17) 

This procedure should be carried out for 0<s<r,l<i< N^^^ and 1 < i < N^'^^ if 

order r is desired. This is an algorithm to obtain the set S{C^^^ * C^'^\ xi^\ xl?^ + A) and has 
been implemented in the python script language, which has a great capability of dealing with 
complicated list operations. In this way, one can obtain the vertices for any parallel transporter. 

3 Extension to the Schrodinger functional 

Let us now proceed to the SF [6j. This section is the central part of the paper. 
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3.1 Preliminary 

We consider, in the following section, a finite box of size xT with periodic boundary conditions 
in the spatial directions. We impose Dirichlet boundary conditions for the link variables at the 
time boundaries, 

U{x,k)\^^^, = eMaC}, U{x,k)\^^^T = exp{aC'}, (18) 

where the boundary fields C, C in the SU(3) gauge theory [Ij are given in Appendix |Bl The 
SF in the theory is defined as 

Z[C,C'] = J D[U]e-^^^l (19) 

It is shown in Ref. [6j that a minimum of the plaquette gauge action is given by the lattice 
background field U{x,ii) = V{x,fi), whose elements are given by 

V{x, 4) = 1, V{x, k) = V{X4) = e*»(^^4-*^), (20) 

where V{x,k) is independent of the spatial direction k. An explicit form of the color electric field 
£ is given in Appendix [Bj For later use, we introduce a convenient notation for the background 
field 

= l-<5^4 = (1,1,1,0). (22) 
By using these notations, we can express the classical link variable as 

v{x,^,r = v{x^r'^^, (23) 

where m takes the values ±1 and V{x,id)"^ represents V{x,^) for m = 1, and V{x,id)^ = 
V{x,fi)~^ for m = —1. 

3.2 An expression for the link variable 

In this section, we introduce a proper expression for the link variable U{x,fi) in the SF case 
which will be useful when generating vertices. In the presence of a background field, the link 
variable is written as 

U{x,fi) = e''3oi>^^'''>V{x,n). (24) 

When dealing with the SF in the presence of the abelian background field, it is convenient to 
use a color decomposition. 



q,{x) = ^qt{x)P, (25) 

a=l 
8 

Qkix) = 5]g^(x)re^<^"(-^)/2, (26) 



a=l 
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Table 1: Phases 4>a{x4)-, 4'a-: 9'q4>a{x4) and ^^(/'a for SU(3) case. The other components a = 2,5, 7, 
which are not shown here, are defined as ^2 = — <^i) (ps = ~04) 07 = ~|?^6- The 7 is given in 
eq. (|133p in Appendix |Bl 



where an overview of the phases 4'a{x4) are given in Table [TJ The idea of inclusion of the phase 
can be found in the original SF paper [6]. The basis of the Lie algebra can be found in [3 [8] 
and is also summarized in Appendix |Al By using in eq. (j22p . the expression for (74 (x) and 
qk{x) are unified as 

8 

q^{x) = ^^^(x)re^'^'"^"(^^)/2_ (27) 

a=l 

By using the above expression, the link variable in eq.(|24p can be written as 



r=l " ai,--- ,ar 

where we have used eq. (j23p for the lattice background field and introduced the notation [II ■ ■ ■ I] 
for representing a 3 x 3 color matrix. 

Now we will derive an expression for U{x,^)'^ = U{x, fi)~^ . Thanks to the phase factor 
introduced in eq. (|28p . it will turn out that it has a form similar to U. Let us show this explicitly. 
The inverse link variable can be expanded in the following way 

U{x,fi)-^ = Vix,n)-^e-''3°i^^'') = V{x4)~''>' 

+ Y T Yl ■ ■ ■ [^(3^4)""^/°' • • • e'"" ^-1 (29) 

This naive form is not suited for our purpose. The location of the V is different from the one in 
eq. (l28p . We want to put V to very right in the color matrix. For this purpose, we employ key 
formulas which will play an important role in the following 

[V{x4)I''V{x4)-^] = [rie*'^'"^^*), (30) 

■a^£ja^-ia^£] = [I''] e''^'- , (31) 

where (p'^ are shown in Table [TJ The formulas are derived from a special property of the Abelian 
background field whose generators are elements of the Cartan sub-algebra. More details on the 



7 



formulas are given in Appendix [Cj Let us return to U ^ . By applying the formula in eq. (|30p to 
eq. (|29p . we obtain the following expression 

U{x,fi)-^ = Vixi)-''^' 

r=l ' ai,---,a,- 

where we have been able to place V in very right in the color matrix. Note the sign of the phase. 
Now, we can express the U and f/^^ in a unified form, for m = ±1, 

Note the location of m. 

In the following the term 'color factor' is reserved for the color matrix part multiplied with 
the phase factors. The corresponding color factor for the single link variable reads 

An essential new ingredient for the case of the SF is this color factor which depends on ^4, m 
and /Li, while a color factor with the zero background field case [/"i/"2 . . . jarj jg independent of 
them. 



3.3 Multiplication of the color factors 

Before entering the case of the SF, it is maybe worth discussing the color factor in the vanishing 
background field case in order to make the difference between the two cases clear. As mentioned 
in Section [2l the color factor in the vanishing background field case is independent of the shape 
of the path, therefore the multiplication of color factors is trivial 

Qai---as Qas+i-'-a,. _ Qai---aaas+i---ar (35) 

On the other hand, for the SF case, due to the presence of the non-vanishing background field, 
the color factor depends on the shape of the path and its multiplication is not as trivial as shown 
above. In this section we will present a multiplication rule for the color factor for the SF. 
Our main finding is that any color factor in the SF of order r can be cast in the form 

C"i-"'-(x«,4, A, B, C, D) = ■ ■ ■ /'^'■y(x/,4)^e*"'^^l eiT.l=i{Cu<f>'a ^ +Duct>a^{xi,^))^ _ (3q) 

_L 'J ^ V ^ 

' 1/(1) phase 



3x3 matrix 



This can be shown by making use of two properties, first, the background field being abelian 
eq. (130|3ip . and second, the linear time dependenc^ of the exponent of V{x4) in eq. pOl) . Note 



In some special case like [9], the property of the linear time dependence is lost. A multiplication algorithm 
for this case is given in Appendix [DJ and it has already been used in a one- loop calculation in the reference. 
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that the color factor has a time dependence which comes from the background field. We choose 
it as x;^4 which is the time component of the coordinate of the left-most link variable of the 
parallel transporter 

P{£,Xs,Xe) = U{xi,firUU-- - , (37) 

with Xs = xi + afi{l — m)/2. In practice, xi is set to the origin. The components A and B are 
single component integer. C and D are multi-component integers of size r 

C = {Ci,---,Cr), (38) 
D = {Di,---,Dr). (39) 

The (p'a in eq.([36]) are the time derivative of the (paixA,) and do them-self not depend on the time. 
Their values are given in Table [TJ In the expression of the color factor, the information about 
the lattice size and the background field parameters (ry, v) is encoded in ^(^^ 4), £, and the 
phases 0a(2;i,4), (p'a- When producing lists, this information is actually irrelevant, only x/^4. A, 
B, C and D are required. The former is only needed when writing down a numerical expression 
of a vertex to a diagram calculation program at a second stage. The benefit of the expression 
is that we can separate the information x;^4, A, B, C and D (which will be included into the 
list structure) and the lattice structure information. Therefore, we can perform a symbolic list 
operation, independently of the detail of the lattice size and the background field parameters. 

Let us give an example of a configuration of {xi^i, A, B,C, D) for the single link U{x,fi)'^ 
whose color factor is given in eq. (|34p 



xi,4 = Xi, (40) 

A = ma^, (41) 

B = 0, (42) 

C = (Or^), (43) 
r terms 

D = {ma^,- ■ ■ ,mafj_). (44) 
r terms 



We have obtained a manageable expression for the color factor. Next we formulate a mul- 
tiplication for the integer list (xi^4^,A,B,C,D). An important point is that even in the SF the 
multiplication of the color factor is closeqj. 

C^l-"^^ (xg , , , C(l) , ) X C^^+i (xg , ^(2) , 5(2) ^ ^-(2) , I)(2)) 

= C'''-''^{xi,4,A,B,C,D). (45) 

From an actual multiplication of the color factors, we find the algorithm to get a correspond- 
ing list from two given lists {x\^Ia(^\b(^\c(^\d(^^) and (xf], A^^), 5(2) , c^^), I)(2)) in the 

^This can be understood from the background field gauge transformation in Ref. [lU] 
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following way 



xi,4: < — x[^l, (46) 

A < — + (47) 

B ^ 5W + 5(2) + At^(2)^ (48) 

C ^ {C[^\- ■ ■ , CP,C[^^ + 25(1) + AtDf\- ■ ■ , + 25(1) + AtDl^},), (49) 

^ V ' ^ ^ ' 

s terms (r - s) terms 



r terms 

D ^ (Dii\...,I)«,Z)f)+2A«,...,Z)(l + 2A«), (50) 

V ' ^ V ' 

s terms (7.-5) terms 



r terms 

where At = x[^) "^Pi • It turns out that the resulting A and B remain single component integers, 
on the other hand, the resulting C (D) is given by combining C^^) (D^^)) and C^^) (D^^)) with 
some additional terms for the latter part. Since x;^4, A, B, C and D are all integer values and 
this operation is simple, the multiplication algorithm is suited for a symbolic operation and is 
easily implemented in the python script language. 

In an actual implementation, as a new ingredient, we have to add the new components xi, 
A, B, C and D to the earlier list structure, 

L('')[/c] = {{ai[k],- ■ ■ ,ar[k]),Xs,Xe,; f[k]} 

{{ai[k],--- ,ar[k]),Xs,Xe,xi,A[k],B[k], 

{Ci[k],--- ,Cr[k]),iDi[k],--- ,Dr[k]y,f[k]}. (51) 

For two given lists, n = 1,2, 

LW[i](£("),x("),x(")) = {(aS"^H,--- ,a(")H),x("),xi'^),xi"\^('^)H,i?(")[i], 

(C;")[i],--- ,C(")[i]),(Z)i")W,--- W);/W[z]}, (52) 

the multiplication algorithm for lists is summarized as follows (the relabelling and amplitude 
parts are the same as before). 

• Relabelling: 

{a^li],--- ,a''p[i\,af\j],--- ,al^l^[j]} — >{ai[k],--- ,as[k],as+i[k],- ■ ■ ,ar[k]}, (53) 

• Color factor part fea.dieD-dSOD): 

(x« , A(i) W , 5(1) [.] , C(i) [.] , Z)(i) H ) * (xg) , ^(2) [j] , 5(2) [,] , C(2) [j] , 5(2) [j] ) 

{xi,,,A[k],B[k],C[k],D[k]), (54) 
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• Amplitude part: 

• Creating the list: 

L(^) H (£(1) , x« , x« ) * L(-«) [j] (£(2) , + A, + A) 

= {{ai[k],--- ,ar[k]),xi^\x^^^ +A,xi,A[k],B[k], 

{Ci[k],--- ,Cr[k]),{Di[k],--- ,Dr[k]y,f[k]}. (56) 

In order to start the multiplication algorithm, we need an initial set of lists. It is the set for 
a single link variable C/(x, /z)™" which is given by the list 

L^"") = {( a,- - ,a ),Xs,Xe,xi,ma^,0, ( 0, • - ,0 ), { ma^, ■ • • ,ma^ );'nf}, (57) 
r terms r terms r terms 

with 

1 — m 

Xs = x + afi — - — , (58) 
1 + m 

Xe = x + afi — - — , (59) 
xi = X. (60) 

To check that the algorithm works properly, we perform a one-loop calculation of the SF 
coupling. The results are shown in Section \5\ By calculating this quantity, we can check the 
two-point vertex. We investigate not only the plaquette gauge action but also the improved 
gauge actions including six-link loops. As a further check we compute the Big Mac diagram, 
which includes three-point vertex, for smaller lattices and check the consistency with the old 
results (private communication with Peter Weisz). 



4 How to exploit the set for the gauge action 

In this section, by using the python output (set 5) generated by the above algorithm, we describe 
how to implement the vertex in a code which calculates Feynman diagrams. We exclusively 
consider the gauge actions, and in both coordinate space and time- momentum space for future 
purposes. In this section we write /x and x explicitly instead of using the super index a = {fi, x), 
and set the lattice unit a = 1. 
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4.1 Coordinate space representation 

The gauge actiorH consists of a loop C 

5G = ^5]Retr[l-t/(C)], (61) 
% c 

where U{C) is the parallel transporter around C. Its expansion in terms of go in the coordinate 
space is given by 

r=2 xi,--- ,Xr 111,--- ,t^r ai,--- ,Clr j=l 

The symmetrized vertex V is given by summation over unsymmetrized vertices Y 

,xr) = ^ a-y;/;:;;(xi,... ,x,), (63) 

where Sr is the permutation group of the order r. The symmetrized vertex may be rewritten as 

'^'-) = ^ E ^■Y;,l}::,:^ixu--- ,xr), (64) 

where the partially symmetrized vertex is 

y'^:.^ (:ei , • • • , X,) = ^ a • y;/.:;-;; (xi , • • • , x,) . (65) 

Zr is defined as a subset of the total permutation group St and consists of the inversion and 
cyclic permutations. The summation over Z^. in eq. (|65p can be done at the python level even for 
the SF. Some details about the partial symmetrization is given in Appendix |El The number of 
elements of Sr and Z^ are = r\ and i^Z^ = 2r respectively. Therefore there are remaining 
permutations ^Sr/Zr = r\/2r = (r — l)!/2. Note that for r > 4 we need to perform the 
remaining symmetrization (summation in eq. (j64p ) at the Feynman diagram calculation stage. 

To make the discussion more concrete, we use the plaquette gauge action as a specific ex- 
ample. First we consider a sum of the plaquette over all directions at fixed point 

^Retr U{0,fi)U{0 + fl,u)U{0 + i',n)''U{0,i^)^ , (66) 

and the corresponding set whose elements are lists of the from 

L^''^[i] = {(wWr-- ,fir[i\),idxi[i\,--- , dxr[i\), 0,0,0, A[i], B[i\, 

(Ci[i],--- ,Cr[i]),m^,--- ,Dr[i]y,m}, (67) 



^We will give more precise definition of the gauge action in Section (5] Here the reality and trace properties 
are important. 
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Figure 1: The Big Mac diagram 



where we assume that the partial symmetrization procedure has been done. As described in 
Ref. [Ij, by making use of the translation invariance in the spatial directions, the dxi, dx^ are 
translated to the standard form [2\, so that the number of lists can be reduced. On the other 
hand, this reduction does not apply for the time direction due to the loss of the time translation 
invariance in the SF. Note that in eq. (|66p the summation over the coordinates is not taken but 
kept. In this case we may set Xs = = xi = 0. If one carries out the sum at the python level, 
it renders the number of lists prohibitively large of order L^T, and this is apparently inefficient. 
Instead the sum is carried out in the diagram calculation code. In order to obtain the complete 
vertex, we have to loop over the xi exhausting the whole elements of the vertex. In this way, all 
non-zero elements of the vertex are assigned. 

Let us show how to use the the set for the partially symmetrized vertex in a diagram calcu- 
lation code. As an example, let us take the Big Mac diagram, which appears in the two-loop SF 
coupling calculation [5] , shown in Figure [H whose expression in the coordinate space is given by 

BigMac = ^ Yl C/.%1(^i'^2,X3)y;/;^^"«(x4,X5,X6) 

{a},{fi},{x} 

X {qZ {xi)qZ {xi)MqZ i^^Kl (^^MqZ i^^Kl (^6)>o, (68) 

where 1/12 is a symmetry factor and V is the symmetrized three-point gluon vertex. The free 
gluon propagator in the coordinate space is expressed in terms of that in the time-momentum 
space D'^j^{p, xo,yo), which gauge fixing is done, 

{q''^ix)qt{y))o = ^^^^.(^^o, yo, x, y), (69) 

and 

-5^1.(2^0, yo,x,y) = ^^^^(xcycx - y) 

= ;^E^'''^""''^^''''''"^'''^^^M-(P'^o,yo). (70) 
p 

The notation a is defined in Appendix |Al In the following, we use the notation t = X4. By 
using the propagator S and performing the partial sum over the color index 04, 05 and ag, the 
BigMac diagram is given by 

BigMac = ^ Yl Yl CA'Ml(*l'*2,t3,Xi,X2,X3)l^™|(i4,t5,i6,X4,X5,X6) 

01,12,03 {fi},{t},{x} 

X S^1m4 XI - X4)5^^^^ (t2, ts, X2 - Xg)^^!^, (ts, te, X3 " Xg). (71) 
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The color factor of order r for the gauge action is given by 



r^^-''^it,A,B,C,D) 



+ {-lyP^ ■■■P'e- 



(72) 



Actuahy almost all elements of the color factor are zero, therefore it is convenient to use a color 
index 

{ai[n],a2[n],a3[n]), (73) 

with running index n which contains indices for the nonzero elements of the color factor. Fur- 
thermore, another list structure is used for the vertex. The symmetrized vertex of order r = 3 
is written as 

!(X1,X2,X3) = lY'%f.f{xi,X2,X3), (74) 



^ 1^1 1^.2 ' 



where Y' is the partially symmetrized vertex. At order r = 3, y and Y' are equivalent up to an 
overall factor 3!. The partially symmetrized vertex can be generated from the python output, 
the list structure for the first gluon vertex in eq. (j7ip is generated from 



(75) 



[i\ = { (^1 [i] , ^2 [i] , W ) , ( {dti [i] , c?xi [i] ),{dt2 [i] , (ix2 [i] ) , (dts [«] , d^s [z] ) ) , 
0,0,0, A[i],B[i],{Ci[i],C2 [i] , Cg [i] ) , (Di [i] , D2 [i] ,D-M);f W } , 



with index i identifying the list. The number of lists is N'i = 396 for the plaquette action. 
For the second vertex we use the index j for identification. The coordinates tm and {m = 
1, 2, 3, 4, 5, 6) are given in terms of the elements of the list as well as t, s, x and y. 



xi = x + dxi [i] 
X2 = x + (ix2 [i] 
X3 = x + dX3[i\ 
X4 = y + dxi [j] 
X5 = y + dX2 [j] 
X6 = y + (iX3 [j] 



h = t + dti[i\, 
t2 = t + dt2\i], 

t3 = t + dt3[i], 
ti = s + dti[j], 

t5 = S + dt2[j], 

te = s + dtslJ]. 



(76) 
(77) 
(78) 
(79) 
(80) 
(81) 



Now the indices (t, x) and (s, y) run over the whole lattice. By making use of the above 
parameterizations and the translation invariance for the spatial direction, the diagram is written 
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as 

^3 52 396 

BigMac = E E EE 

72(00101 loop) ^,j(python loop) t,S X 

X '"''^'"''^'"^ ^b1, 5b1, 

X wli bl (i + 1 W , 5 + 1 [j] , X + dxi [i] - dxi [j] ) 

^ wL [J] it + dt3[^,s + dts [j] , X + dx3 [i] - dx3 [j] ) . (82) 

Note that the summations over {fi}, {t} and {x} in eq. (|68]) of order 4*^ x T^{L^f are replaced 
by i,j (python loop), t,s and x of order (396)^ x T^Li^. In the summation step, one has to be 
cautious about the boundary conditions, 

qk{x)\xi=Q = qk{x)\xi=T = (83) 

Here we give a short comment about the color factor of the gauge action. For any gauge 
action composed of closed and traced loops, it holds that A = Q. This means that the color 
matrix part (which is the color factor without the phase factor that depends on X4), 

tr [/"I • • • pr^iSB ^ (_l)r-ja. . . . jaig-i^Bj ^ (g4) 

is time independent. Therefore one does not have to calculate the color matrix for every X4 when 
writing down the vertex in computer code at the Feynman diagram calculation stage. Before 
the time loop X4, one can store the color matrix in somewhere, and then in the time loop it can 
be called. This is realized by virtue of separating the color factor into two parts, the information 
about lists and the background field. 



4.2 Time-momentum space representation 

The Fourier transformation for the quantum field is given by 

= 73 E ^^^■^"^^''''^^m(P' ^4), (85) 



L3 
p 



where the spatial momenta are given by pk = ^iruk/L with = 0, L — 1. An expansion of 
the action in terms of the coupling go in the time-momentum space is given by 

r=2 ■ ^ ^ Pi,---,Pr 

r 

12 V^-::^:;{pu--- ,Pr-,tu--- ,tr)YlqZ{-Pj,t,). (86) 

/il,--- ,/ir ^l,-"" 0.1, ■■■ ,ar J — 1 
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Similar way to the case in the coordinate space, the symmetrized vertex is given by the 
partially symmetrized vertex 

^■■.■;;(Pl'--- 'P-*l'--- = ^ E ^•^;i^.:^r(Pl,--- ,Pr;tl,--- ,tr-). (87) 

We also have written python code for the time-momentum representation. The list structure 
for the time-momentum space is a little bit different from the one in the coordinate space. The 
differences occur in the coordinate of the link variables only, x — > (t, v), where i = X4 and v is 
a spatial component of mid-point of the link variables in two lattice units. The list for eq. (|66p 
is given by 

i^'^W = {(mWr-- ,/i.W),((tiW,vi[i]),--- ,(t,,[i],v,[i]), 0,0,0, 

yl[i],SH,(Ci[i],--- (88) 

By making use of the propagator in the time-momentum space 

kfykv, 3;o)^!!(q, yo))o = '56a^^'5p+q,o-C>|^;.(p, a;o, yo), (89) 
the Big Mac diagram in the space is written as 

^ 52 396 

BigMac = E E EE 

n(color loop) ij' (python loop) ^i* Pi iP2 

X l/[j]C^iN^2N^3N(,^ C[j], Z)[j])el(P--^b-]+P--^yi-^-3b1) 

xDlf^^^^.^iq, t + dts[i],s + dtsbD, (90) 



where qk = 2-irmk/L with = nik + n2k mod L for pik = 2-Knik/L {i = 1,2). We have used 
the relation, 

D%{V, xo, yo) = D^^i-p, yo,xo). (91) 

Note that the computational effort scales with T^L^ which is more demanding compared to the 
coordinate space by L^. Even though the cost for the Fourier transformation of the propagator 
in eq. (j70p scales proportional to T^L^, the prefactor is negligible compared to the sum in the 
diagram whose factor is 52 x (396)^. This shows the benefit of using the coordinate space for 
the computation of bubble-type diagrams [5]. 
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5 Applications 



In this section, we show some appHcations of the algorithm, especially we consider the one-loop 
computation of the SF coupling. From this computation we can check that the algorithm works 
properly and furthermore we can investigate more complicated gauge actions, like improved 
gauge actions including chair type and parallelogram loops. 



5.1 Improved gauge actions 

In this subsection, we summarize the definition of the improved gauge actions including six-link 
loops. More details about the definition can be found in [11]. The improved gauge action for 
the SF is given by 

2 ^ 

'^G = - E E Wi{C,gl)ReiT [1 - U{C)] , (92) 
^0 i=oce5i 



where U{C) is a parallel transporter along a loop C. Since we are taking the real part, loops 
C which differ by orientation only are considered equal. Si denotes sets of loops C on the 
lattice as given in Figure [2] together with the corresponding factors Cj , which are normalized by 
Co -I- 8ci -I- 16c2 + 8c3 = 1. We adopt the weight factors for "Choice B" in Ref. [II]. For the 
plaquette, this means 



and for rectangular loops 





if C lie on one of the boundaries, 

if C just touch one of the boundaries, 

otherwise. 



if C lie completely on one of the boundaries, 
if C have exactly 2 links on a boundaries, 
otherwise. 



(93) 



with the perturbative expansion of the boundary counter term [T2] 

Cc^gl) = co{l + c^^'\l + 0{gt)\ 
c,c^{gl) = c,{2,/2 + c^^'\l + 0{gt)). 



(94) 



(95) 
(96) 



The tree values of the boundary counter terms were given in Ref. [TI]. The one-loop coefficients 
for the improved gauge actions composed from plaquette loops and rectangular loops (Iwasaki, 
Symanzik and DBW^j) were already determined in Ref. [T2] with a constraint 



2c, 



(97) 



As discussed in the reference, this constraint is convenient when we deal with one-loop 0(a) 
improvement, and we assume this in the following. As we shall see later, the remaining one- 
loop values c^ will be fixed for the other improved gauge actions including the chair and 



'Note that the value of ci for the DBW2 here is a bit different from that in |12| . 
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Co Cl C2 C3 

Figure 2: Loops together with corresponding factors. 



parallelogram loops. The weight factors for the chair and parallelogram loops (i = 2, 3) are 
given by 

,^ 2\ f if C lie completely on one of the boundaries, , . 

W,{C,go)-^ c, otherwise. ^^^^ 

Gauge actions including the chair and parallelogram loops are hardly used in current simula- 
tions, but here we consider such actions for comparison with the old perturbation theory results. 
We deal with seven gauge actions whose name and values of Cj are given in Table [2j 

5.2 One-loop results for the SF coupling 

We perform the one-loop computation of the SF coupling in pure SU(3) gauge theory for the 
improved gauge actions. We set T = L as usual. The SF coupling is defined [1] by 



(99) 



r)=u=0 



with r = — In ^ [C, C] and a normalization 

k = l2[L/af [co{sin(27) + sin(7)} + 4(ci + 2c2 + C3){sin(47) + sin(27)}] , (100) 

where 7 is given in eq. (|133p . 

The coupling constant has the perturbative expansion 

-glAL) = gl + m^l\L/a)gt, + 0{gl), (101) 

and the one-loop coefficient is given by 

m^^\L/a) = mf\L/a)-2c^'^^^R{L/a)a/L, (102) 
^ ^ co{sin(27) + sin(7)} + 4ci{sin(47) + sin(27)} 

^ co{sin(27) + sin(7)} + 4(ci + 2c2 + C3){sin(47) + sin(27)} ' 

mj^\L/a) is the one-loop coefficient which is computed with the tree level value of the boundary 
counter terms, that is, the tree level 0(a) improved coefficient. On the other hand, m^\L/a) 
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Action 


Cl 


C2 


C3 


^0 


Ai 


Plaquette 











0.3682819(6) 


-0.1779(3) 


Wilson RG 


-0.252 





-0.170 


-0.2184103(8) 


0.2987(3) 


Iwasaki 


-0.331 








-0.204903(3) 


0.304(2) 


DBW2 


-1.4088 








-0.62816(5) 


0.90(2) 


Symanzik 


-1/12 








0.1361508(3) 


-0.0060(1) 


Symanzik II 


-1/12 


1/16 


-1/16 


0.1403119(5) 


0.0216(2) 


Symanzik III 


-1/12 


-0.1 


0.1 


0.1491183(6) 


-0.0745(3) 



Table 2: The coefficients and Ai of asymptotic expansion of the one-loop coefficient for the 
various gauge actions. 



is a one-loop 0(a) improved one which is expected to have no a/L term. What we compute 
actually is mf'\L/a) in a range L/a = 6, ...,64. The contribution is separated into two parts, 
the gauge field and the ghost parts. 



(0) 

2k 

a,p 



-^VTr 



+ (ghost contribution). (104) 



The inverse propagator i^"(p) is a sum of the two-point vertex V^l^{p, —p;t, s) in eq. (|86p and 
the conventional gauge fixing term [3]. We use the ghost action in the previous reference. The 
rj derivative of K'^{p) can be build numerically from the python output as shown in Appendix 
[Gl In the one-loop computation we adopt the time-momentum representation which is still 
advantageous at this order. The trace is taken over the Lorentz and time indices. Cp"* in 
eq. (jl02p can be chosen to cancel the linear term of a/L in 'mf'\L/a), and then m^^\L/a) turns 
out to be the 0(a) improved quantity. 

By analyzing the tree level 0(a) improved one-loop results following [13], we obtain the first 
several coefficients of the Symanzik expansion form. 



,(0) 



(L/a) = Ao + Bo ln(a/L) + Aia/L + Bia/L ln(a/L) + 0{{a/Lf). (105) 



Bq should be 26o, which is the universal one-loop coefficient of the beta function. We check that 
|i?o/(2&o) ~ 1| < 10~^ holds for all improved gauge actions, which we calculate. We also observe 
that \Bi\ < IQ-^ for all cases, and i?i = is a signal for the achievement of the tree level 0(a) 
improvement. By assuming that Bq = 26o and i?i = we obtain the values of and Ai, which 
are given in Table [2j Their values for the case of the plaquette action (the explicit value 
of Aq is shown in [iSj), and rectangular type actions [12j . (Iwasaki and Symanzik) are already 
known and we confirmed the values to ensure the consistency. 

Let us discuss the ratio of the Lambda parameters. Actually there are several old pertur- 
bation theory results of the ratio of the lambda parameters between the plaquette gauge action 
and the improved gauge actions, therefore we can compare with our results. Before taking a 



19 



continuum limit, we have to perform the renormalization. We introduce a renormalized couphng 
5Lat(/^) through 

al = 9Lt{f^) + zi{a^)-giM + 0{9iM), (106) 
zi{a^i) = 2&oln(a/i). (107) 

By substituting eq. (jl06p into the expression of g^^{L) of eq. (jlOip with eq. (jl05p . we obtain 

fsAL) = -glM + -glM [Ao + i?o ln(a/L) + 2ho ln(a/i) + 0{a/L)] + 0{glM)- (108) 

By noting Bq = 2bQ and setting /x = 1/L, we can cancel the log divergence. Then we can take 
the continuum limit 

ghiL) = gLdVL) + AogtdVL) + 0{gl,,il/L)). (109) 
The ratio of the lambda parameters is given in terms of Aq 

ALat 



AsF 

From this, we can compose a ratio 



exp[-Ao/(26o)]. (110) 



Ai^ = ^ilBP£^ = eMiAl'^"" - Ar')/i^bo)]. (Ill) 

Apiaq Apiaq/ AsF 

The resulting values are given in Table O we observe rough consistency with the old values. 
Another ratio of the lambda parameters in the "lattice scheme" and the MS scheme is given by 

ALa. ^ ^^ = ^^p[_^A, + c^/{47r))/{2bo)], (112) 

^MS ^SF A^jg 

ci\Nt=o = 1.255621(2), (113) 



where ci is taken from the two- loop SF coupling paper |T3j , its value for the several gauge actions 

is given in Table [H This table may be useful for future references. 

From the resulting value of Ai we can determine the one- loop coefficient of the 0(a) boundary 
p(i) 

counter term . The improvement condition yields 



p(i) _ Ai 



(114) 



' 2(co + 8ci)' 
Together with eq. (|97p . one can achieve one-loop 0(a) improvement. 

5.3 Relative deviation of the step scaling function 

Finally, we investigate the relative deviation of the step scaling function (SSF) for the various 
gauge actions. The SSF for the running coupling [21] is defined by 

a{u) = glp{2L), u = 5sf(^)- (115) 
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Action(impr.) 


our results 


[14J 


[15J 


m 


m 


m 


m 


Wilson RG 


67.4384(5) 






67.97(9) 




67.6(3) 


67.37(14) 


Iwasaki 


61.207(2) 


61.2064 






59.05(1.00) 






DBW2 


1277.1(5) 


1276.44 












Symanzik 


5.29209(3) 


5.29210 


5.294(4) 


5.29(1) 




5.29(1) 


5.29(1) 


Symanzik II 


5.13636(3) 










5.13(1) 




Symanzik III 


4.82173(3) 






4.842(2) 









Table 3: Status of Aimpr/Apiaq for pure SU(3) gauge theory. The results from Ref. [T6l[T7] tend 
to disagree with ours. Other values from [TU [TSJ HSl [12] are consistent with our results. Note 
that the value of ci in [T3] is a little bit different from ours and actually they use ci = —1.4086. 



Action 


^Lat/A^jg 


Plaquette 


0.0347109675049892(2) 


Wilson RG 


2.34086(2) 


Iwasaki 


2.12455(5) 


DBW2 


44.33(2) 


Symanzik 


0.1836938(4) 


Symanzik II 


0.1782883(4) 


Symanzik III 


0.1673674(8) 



Table 4: ALat/Aj^ of SU(3) gauge theory. For the plaquette gauge action it is taken from 
Ref. [2D], where a direct calculation of the ratio was done. Our results are those for the other 
gauge actions. 



We denote the SSF measured on the lattice with Y!,{u,a/L). It is expected to converge to the 
continuum value, 

a(u) = lim S(n,a/L). (116) 
A scaling behavior to the limit can be described by the relative deviation which is defined from 

d{u,aL) = . (117) 

a{u) 

In perturbation theory, it is expanded as 

6{u,a/L) =6i{a/L)u + 0{u'^). (118) 
The one-loop coefficient is given by 

6[^\a/L) = {2L) - mf\L/ a) - 2boln2, (119) 
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(1) 
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Wilson RG ■— 

Iwasaki 

DBW2 

Symanzik — 
Symanzik II 

Symanzik III 
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-0.03 
-0.04 
-0.05 
-0.06 
-0.07 



0.005 0.01 0.015 0.02 0.025 0.03 
(a/L)' 



Figure 3: Left panel is for the tree level 0(a) improved case, and right one is for the one-loop 
relative deviation where the 0(a) term is removed. 



where k = 0,1 denotes the degree of improvement, the results are shown in Figure [3] for the 
various gauge actions. Symanzik type actions have relatively good scaling behavior, while the 
renormalization group improved types show rather large cutoff effects at one-loop level. 

6 Concluding remarks 

In this paper we derive the multiplication algorithm to generate Feynman rules for the SF with 
an abelian background field. The essential new ingredient for the extension is how to treat the 
color factor which involves the background field. By making use of the key formulas eq. ijSOp and 
eq. ([3T|) . which are derived from special properties of the abelian background field (eq. (jl37p ). this 
difficulty is resolved. We present the multiplication algorithm for the color factors in eq. ()46t - 
[50]) . By making use of the algorithm one can obtain vertices to any order for any shape of the 
parallel transporter in the SF. One has to keep in mind that as shown in Appendix|E]one needs 
to perform the remaining symmetrization for the vertices r < 4 of the gauge action in order to 
obtain the totally symmetrized vertex at the Feynman diagram calculation stage. 

As a first application and to check the correctness of the python code, we calculated the 
one-loop coefficient of the SF coupling for various gauge actions. We observe consistency by 
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comparing with previous calculations, and analyze the scaling behavior of the step scaling func- 
tions. As a further check we compute the Big Mac diagram for smaller lattices and observe a 
consistency with the old results. As a further application, we applied the automatic method to 
a lattice with L = T ziz a which is motivated by considering staggered fermions in the SF [9] . 
Actually, on that lattice, the time dependence of the background field is not uniform anymore. 
Therefore we extend our algorithm to this case, the details of the extended algorithm are shown 
in Appendix [DJ A preliminary result for the gauge sector by making use of the algorithm was 
already presented in [9^ i22j . 

It is desirable to go beyond one-loop for actual applications. In that line, the full two-loop 
computation of the SF coupling is the next target. In this paper, we exclusively consider the 
gluonic sector as a specific example. Of course, however the multiplication algorithm can be 
applied for the fermion actions, for example, Wilson fermions and also for the clover term. The 
algorithm is suited for all link connected objects. Therefore an application for the HQET with 
highly smearing is also possible. 

Some readers, who are interested to use the generalized version of the python script, are 
invited to ask the author to send it. 
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A The basis of the Lie algebra su{3) 



For the Lie algebra su{3) one may choose as a basis 

with the Gell-Mann matrices A'^ for a = 1,2,..., 8, such that 



tr 



1 



•ab 



rpa rpb 



ahcrjic 



jaOCrpt 



(120) 



(121) 



For the SF case, however, it is convenient to choose another basis. New matrices A" [7l|8] are 
introduced, which coincide with the Gell-Mann matrices A" except for the two diagonal matrices. 



A3 
A^ 



-^A^ + ^A^. 
2 2 

2 2 



(122) 
(123) 
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The matrices A" are 



A^ 



A^ 



/ n 


1 

1 


n 


u 


U 


U 















1 


























(: 





1 




1 













1 












-1 





—i 





u 


U 


U 















—i 


























(: 





—i 




i 







' 2 





71 ( 



. 


-1 





After a normalization, 



rpa \ a 



these matrices may be used to define a new basis which is given by 









/2 


= r_ 














= U- 






/6 






f 


= V- 







for the non-diagonal matrices and 

for the diagonal ones. For this basis, one has 



(124) 



(125) 



(126) 



(127) 



(128) 



where 1 = 2, 4 = 5, 6 = 7, and 2 = 1 and so on. For the diagonal matrices, one has 3 = 3 and 
8 = 8. The normahzation is chosen such that 



tr 



jajb 



2 



(129) 
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B The boundary fields and the constant color electric back- 
ground field 

In the basis chosen in Appendix [Aj the boundary fields are expressed as 



C 



I 

T 

i 

L 

i 

L 

i 

I 



7? - f 

r?(-i+zv) 



-7?(i+zy) + f 

-r? - vr 

+ ^ 

r]{\-u) + ^ 



(r? + 7r)A^ 



The constant color electric background field £ is proportional to 



2 
£: = -7 I -1 
0-1 



with 



1 / vr 



The ?7 derivative of £ is given by 



^ = -2./3.^/« = ^ 
or] or] LT 



(130) 



(131) 



(132) 



(133) 



(134) 



C The Cartan sub-algebra 

Let us assume that Hi for i = 1, - ■ ■ ,m [m is the rank of the algebra) are hermitian and elements 
of the Cartan sub-algebra. They commute, 



The generators Ea of the original algebra satisfy 

[Hi, Ea] = ]J-aiEa, 



(135) 



(136) 



where ]j,ai are the roots, which are the weights of the adjoint representation. By making use of 
eq. (ll36p . one can see that 



(137) 
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a 


1 


2 


3 


4 


5 


6 


7 


8 




i 

2 


i 

2 





i 

2 


i 

2 


1 


-1 





fJ-a2 


2 


V3 
2 





2 


2 












Table 5: The roots fiai for SU(3) case. 



with real coefficients hi. An Ea sandwiched between group elements generated by the Cartan 
algebra does not mix with other Eh {b ^ a), but turns out to be the Ea itself multiplied a phase 
factor. 

In our case of su{3), the Cartan generators are identified as Hi = il"^, H2 = il^, and the 
other generators as Ea = The values of the roots for the sn(3) case are shown in Table [5l 
In terms of the Cartan generator, the background field in eq. (|20p can be written as 

1/(2:4) = e*S^■^^(^^)■^^ (138) 

where the coefficients hj^x^) can be read off from the definition of V^x^). By making use of the 
above formula, one can see that the phase (paix^) in eq. (j30p is given in terms of the coefficients 
hj and the roots fiaj 

(paiXi) = '^hj{xi)lJ.aj. (139) 
j 



D Multiplication for a background field with an arbitrary time 
dependence 

In this section, we provide a multiplication algorithm for the color factor for an abelian back- 
ground field which has arbitrary time dependence. We start from a general form of the color 
factor of order r for this case, 

V''nim) '■max) ^) -'-^y 
= [la, ■ ■ ■ /a.^(tmin)^'-'"^(tmm + l)^*-i„+i . . . y(W)^'-ax] x 

X e3 Eu = l (-^".tmin'^''" (tmin) + -Du,t^;^ + l<^a„ (tmin + l)H \-Du,tni!,^<t>au (tmax)) _ (140) 

For a certain configuration {ai, • • • , 0,-}, the color factor is specified by the set {tmim ^maxj ^) D}. 
tmin and tmax are defined in Figure HI The explicit forms of A and D are given by 

A = (A 

D = (A„,,„A_+i,--- ,A._), (142) 

A = {Dlt, D2t, • • • , At), for tmin < t < tmax. (143) 
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t 



^max 
^end 
^start 



Figure 4: A line which show t^ax, imin, istart and iend- 

Note that all elements of A, D are integer valued. Let us sec an example of (A, D) for the color 
factor for the single link variable U (x, /u)™. For a; = and /x = 4 (tmin = 0, tmax = 1), 

{A,D) = {{Ao,A^),{Do,D^)) 

= ((0,0),( (0,0,^-,0)^ , (0,0,^...,0); ). (144) 

r terms r terms 

For a; = and /x = 1, 2, 3 (tmin = 0, tmax = 0), 

iA,D) = {Ao,Do) 

= (0,(0,0,...,0)). (145) 
r terms 

Even for this case, the color factors are closed under the multiplication (r = ri + r2), 
C"-"'^^ (t!^L, t'^, ^(1), Z)(l))C»n+i-n+^. (t(^)^, t(^l^ d(2)) 
= C»i-"'-(tn,i„,W,^,i^)- (146) 
So we found the following multiplication algorithm. First, the resulting tmin and imax are given 

by, 

t^in = min(t2L,iSn), (147) 
W = max(t(^L>«^^L)- (148) 
Before operating A and D, we have to pad zeros, from the left 
^(i)^i(i) = (if) ) 

^ 'rnin t-max ■' 



( 0, 0, , A% , ), for > tmin, 



(1) 



.(1) 



(^min ~ *min) tCrmS 



(149) 



(AC-) 4(1) \ 

l^.(l) , •■•5^,(1) j, 



for f^^^ — t ■ 

^•-■^ ''min ~ "-mm, 
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where we note that the size of A^^^ is tmax — imin + 1, while that of yl*-^-* is tmL — tmin + 1- Also 
from the right, 

^(i)^i(i) = (id) ) 

^ t-min '•max 



''max 



(tmax - tmL) terms (150) 
/l(^) ^ forf(^) -t 

^ 4-W ' ,(1) /' '■max — ■'max- 



In a similar way, we can pad the others with zeros and obtain A^'^^ and D^^'"^^ . Then the operation 
is done for the tilde objects, 

At = ir^+if\ (151) 

A = (5ii\---,5i;l,5i? + 2i«,...,^gJ + 2if)), (152) 



ri terms r2 terms 



r terms 

for imin < t < tmax- For the gauge action, the symmetrization and the reduction of the number 
of lists can also be applied in this case in a similar way to Appendix [El and iFl 

E Partial symmetrization of the vertex at the python level 

For the gauge action (traced and real), the summation over all possible permutation in eg. (1870 
can be done partially at the python level, that is, by creating new lists corresponding to a 
symmetrized vertex. It is maybe worth starting from the vanishing background field [HE]. The 
color factor in this case {£ = 0, (pa = (p'a — 0) is given by 

c^i-'^'- = trii"^ ---r- + {-lyp- . . . ri]. (153) 

This is an eigenstate of the inversion and cyclic permutations. 

a ■ Cg = Xr{cT)CG, (154) 

where the eigenvalues are given by 

, , f 1, for cyclic permutations, /-^..n 

Xrifj) = < , ,\r f ■ , 155 

[ (~1) 7 for inversion permutation. 
By using these properties, the symmetrized vertex can be written a^ 

V^,\-::^:{pi,--- ,Pr;ti,--- ,tr) = ^_ a-C^i-'^'-a-y^...^^(pi,--- ,p,;ti,--- ,i,), (156) 



^ Here we show the time-momentum space representation, but the following discussion is also valid for the 
coordinate space case. 
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with 



with 

^A-M,(Pl> • • • ,Pr;tl,--- ,U) = ^ Xr(o-)o--y^^...^^(pi,--- ,pr;h,--- ,tr). (157) 

In Ref. [2], it is described how to perform the summation in eq. (jl57p at the python level. We 
have seen the nice properties of the color factor of the gauge action in the case of a vanishing 
background field. At first sight, however, it is nontrivial that the color factor for SF has such 
properties. We will discuss this issue in the following subsections, and explain how to obtain a 
set for the partially symmetrized vertex. 

E.l The inversion permutation for the color factor 

First, let us see the inversion permutation for the color factor of the gauge action in the SF. An 
inverted color factor is written by (remember that A = for the gauge action) 

fTinv-C''""''""'''(^4,0,5,C,L») = C^'^'--'-''^'''{x4,0,B,C,D) 

= {-IYC^'"'-''''{X4,0,B,C,D), (158) 

B = -B, (159) 

(Ci, C2, • • • , Cr._i, Cr) = (Cr, Cr-1, • • • , C2, Ci), (160) 
{Di,D2,--- ,Dr-l,Dr) = {Dr,Dr-l,--- ,D2,Di). (161) 

Under the inversion permutation, the color factor is not an eigenstate anymore. However the 
structure is still stable but with different components B, C and D. 

E.2 Cyclic permutations for the color factor 

Secondly, let us consider the cyclic permutation. A cyclically permuted color factor is written 
as 

acyc-Cir'""''"'''''{x4,0,B,C,D) = C^2'^«-'^'-'^i(x4,0,S,C,L>) 

= C^'"'-"'''{x4,0,B,C,D), (162) 

where we have used eq. ([3T]) . The explicit forms of the "hat" objects are as follows, 

B = B, (163) 

{Cl,C2,--- ,Cr-l,Cr) = {Cr-2B,Cl,C2,--- ,Cr~l), (164) 

{Di,D2,--- ,Dr-i,Dr) = {Dr,Di,D2,--- ,Dr-l). (165) 

Composite cyclic permutations can be done in a similar way. For 1 < j < r — 1, 

{acycY ■C^'-''^{x4,0,B,C,D) = C^+'""''-''''-'''''''"^'{x4,0,B,C,D) 

= C^'"'-'^'-{x^,0,B,C,D), (166) 
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with 

B = B, (167) 

(Ci, C2, • • • , Cr) = {Cr-j+1 — 2B, Cr—j+2 — 2i?, • • • , Cr — 2B, Ci, C2, • • • , Cr-j), (168) 
{Dl,D2,--- ,Dr) = {Dr-j+l,Dr-j+2,--- ,Dr,Di,D2,--- ,Dr-j). (169) 

As a further example, we show an opposite cychc permutation. 

{acyc)~^ ■ C^'"'--"''- {X4,0, B,C, D) = C^''''-''^-"''-'{x4,0,B,C,D) 

= ^^"''■■■"^-'"'■(2;4,0,^,C',D), (170) 

with 



B = B, (171) 
(Ci, (72, • • • , Cr-i, a) = {C2, C3, • • • , Cr, Ci + 2i?), (172) 
{Di,D2,--- ,Dr-l,Dr) = {D2,D3,--- ,Dr,Di). (173) 

If one uses eq. ()194p . it is not so difficult to see that 

(acyc)-' ■ Cg = {acycY'^ ■ Cg. (174) 
E.3 The partially symmetrized vertex for SF 

Let us see how to construct the set for the partially symmetrized vertex for the SF, 

C---m""(pi' ■■■ ,Pr;h,-- - ,tr)= Y^a- y;,\::-;;(pi, • • • , p.; h,--- , u). (175) 

Starting from lists for the unsymmetrized vertex L('')[i], one can create corresponding lists in 
the following way, 

^yf^ (cT,y,)i.LM[i] (l<i<r-l) 

inverse i.^'^^) 

dinv • L(^) [i] {a,y,y . dinv • L(''-) H (1 < j < r - 1). 

Note that in the SF case the cyclic and inversion permutations are not commutative. An 
algorithm to obtain an inverted list, cJinv ■ ^^'^^[^l, from the original list L('')[i] 

LMf,] (177) 
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is as follows, 

{f^l,fJ-2,- ■ ■ ,IJ-r-l,fJ'r) {fJ'r,fJ-r-l,- ■ ■ ,fJ-2,fJ'l), (178) 

{ti,t2,- ■ ■ ,tr-l,tr) > (ir, • • • 1 *2, il), (179) 

(Vi, V2, • • • , Vr_l, Vr) > (Vr, Vr_l, • • • , V2, Vi), (180) 

B — > -B, (181) 

{Cl,C2,- ■ ■ ,Cr-l,Cr) > {Cr-jCr-l,- ■ ■ ,C2,Cl), (182) 

{Di,D2,--- ,Dr.i,Dr) {Dr,Dr-i,--- ,D2,Di), (183) 

/ (-1)7- (184) 

For a cyclically permuted list a eye ■ L^'''^[i], 

^(^•^[i] ^fTcyc-LW[i], (185) 

it is given by, 

(Aii,/i2,--- ,/"r-l,/"r) ' (/X2,/U3, • • • ,;U,.,/Xi), (186) 

{ti,t2, ■ ■ ■ ,tr-l,tr) > {t2,t3,- ■ ■ ,U,tl), (187) 

(Vi, V2, • • • , Vr-l, Vr) > (V2, V3, • • • , Vr, Vi), (188) 

B — > B, (189) 

(Ci, C2, • • • , Cr-l, Cr) — > (C2, C3, • • • , Cr, Ci + 2i?), (190) 

(Z)l,I)2,--- ,^r-l,I)r) (Z)2,I)3,--- ,^r,A), (191) 

f f- (192) 

By gathering all lists, [i]}, {^inv • i^^"^) W}, {(^cyc)^' • L^^''^ [i] 1 1 < j < r - 1} and {{aejcV ■ tiinv • 
L('')[z]|1 < j < r — 1}, the set S" for the partially symmetrized vertex can be composed. 



F Reduction of the number of lists 

One can reduce the number of lists by making use of the properties of the background field. 
From eq. (j30p and (I31h . we found that the color factor obeys the relation 

C^'-''^{xi,A,B,{Ci,--- ,Cr),{Du--- ,Dr)) 

= Q-^^{x4,A,B,{Ci+a,--- ,Cr + a),{Di+p,--- ,Dr + f3)), (194) 

for arbitrary real numbers a, (3. This equation tells us that one can shift all elements of the 
C = (Ci, • • • , Cr) and D = {Di, ■ ■ ■ , D^) to a standard form such that Ci ^ and Di — > 0, 
that is, choosing a = —Ci and (3 = —Di. After performing the shift for C and D of all lists 



31 



in a set, there might be some hsts (e.g. two hsts ^(''•'[i] and -^^^^^[j] with i ^ j) which have an 
identical configuration C[i] = C[j] and D[i] = D[j] (of course the Lorentz index configuration, 
the coordinate configuratiorj^, A and B in the two hsts should be identical beforehand, but their 
amplitudes do not have to be so). In this case, we can merge the two lists into a list by summing 
the amplitude factor f[i] + f[j]. In this way one can reduce the number of lists. Actually this 
achieves a few percents reduction of the size of a set for typical gauge actions. 



G T] derivative of the vertex 

When one calculates the SF coupling, the rj derivatives of the vertices are required. We can 
easily obtain them from the same python output for the usual (non ry-derivative) vertices. Since 
the information of the background field {V{xi), £, cp'^ and cpaix^) which have rj dependence) 
is encoded in the color factor, we need to consider the color factor only for our purpose. For 
example, the explicit form of the rj derivative of the color factor for the gauge action (A = 0) is 
given by 

dQ-''^ix4,0,B,C,D) 
dr] 



LT 



tr 



rai 



I"- re' 



r t8 ja,. 



-lypT 



^EB 
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+ lY. + I)«5,0,„(X4)) C^i-'^^(X4, 0, B, C, D), 



(195) 



u=l 



where we have used eq. ()134p . The information oi B, C and D is enough to identify the above 
term. When writing down the rj derivative vertex on a diagram calculation code, one needs 
additional phases dri4>'a d'q4'a{xi) which are given in Table [1] in analytic forms, and an 
additional color matrix 



tr 



jai _ _ _ jar jSi^ia?€B ^ \y I""^ ■ ■ ■ /"^ig"**^^^^ 



(196) 



This has a form similar to that in eq. ()84p of the usual vertex apart from the presence of I® and 
the relative sign between the two terms. 



References 

[1] M. Liischer and P. Weisz, Efficient numerical techniques for perturbative lattice gauge 
theory computations, Nucl. Phys. B266 (1986) 309. 

[2] A. Hart, G. M. von Hippel, R. R. Morgan and L. C. Storoni, Automatically generating 
Feynman rules for improved lattice field theories, J. Comput. Phys. 209 (2005) 340-353 
[hep^lat704lT026| . 

* In the momentum space representation, the replacement x[i] {t[i],v[i]) is understood. 



32 



[3] A. Hart, G. M. von Hippel and R. R. Morgan, Leptonic widths of heavy quarkonia: S-wave 
QCD / NRQCD matching coefficients for the electromagnetic vector annihilation current 
at 0(alpha(s) v**2), jhep-lat/0605007[ 

[4] M. Liischer, R. Sommer, P. Weisz and U. Wolff, A precise determination of the running 
coupling in the SU(3) Yang-Mills theory, Nucl. Phys. B413 (1994) 481-502 
[hep-lat/93090 05|. 

[5] R. Narayanan and U. Wolff, Two loop computation of a running coupling lattice Yang- 
Mills theory, Nucl. Phys. B444 (1995) 425-446 | hep-lat/950202lj . 

[6] M. Liischer, R. Narayanan, P. Weisz and W. Ulli, The Schrodinger functional: A 
renormalizable probe for nonabelian gauge theories, Nucl. Phys. B384 (1992) 168-228 
|hep-lat/9207009|. 

[7] P. Weisz, "csw note (private note)." 1996. 

[8] S. Kurth, "The renormalized quark mass in the Schrodinger functional of lattice QCD: A 
one loop calculation with a nonvanishing background field." 2002. 

[9] P. Perez-Rubio and S. Sint, The SF running coupling with four flavours of staggered 
quarks, Po^ LATTICE2007 (2006) 249 [0710 . 0583]. 

[10] U. Wolff, "Symbolic tools for perturbation theory for the Schrodinger functional." 2007. 

[11] S. Aoki, R. Frezzotti and P. Weisz, Computation of the improvement coefficient csw to 

1-loop with improved gluon actions, Nucl. Phys. B540 (1999) 501-519 [hep-lat/ 9808007|. 

[12] S. Takeda, S. Aoki and K. Ide, A perturbative determination of 0(a) boundary 

improvement coefficients for the Schroedinger functional coupling at 1-loop with improved 
gauge actions, Phys. Rev. D68 (2003) 014505 [hep-lat/0304013] . 

[13] ALPHA Collaboration, A. Bode, P. Weisz and U. Wolff, Two loop computation of the 
Schroedinger functional in lattice QCD, Nucl. Phys. B576 (2000) 517-539 
[hep-lat/9911018|. 

[14] A. Skouroupathis and H. Panagopoulos, Lambda-parameter of lattice QCD with Symanzik 
improved gluon actions, Phys. Rev. D76 (2007) 114514 [0709.3239J. 

[15] P. Weisz and R. Wohlert, Continuum limit improved lattice action for pure Yang-Mills 
theory. 2, Nucl. Phys. B236 (1984) 397. 

[16] Y. Iwasaki and S. Sakai, The lambda parameter for improved lattice gauge actions, Nucl. 
Phys. B248 (1984) 441. 

[17] Y. Iwasaki and T. Yoshie, Renormalization group improved action for SU(3) lattice gauge 
theory and the string tension, Phys. Lett. B143 (1984) 449. 



33 



[18] A. Ukawa and S.-K. Yang, Asymptotic freedom scales for SU(N) lattice gauge theory with 
improved action, Phys. Lett. B137 (1984) 201. 

[19] W. Bernreuther, W. Wetzel and R. Wohlert, Lambda parameters for lattice Yang-Mills 
actions containing plaquettes with six links, Phys. Lett. B142 (1984) 407-409. 

[20] M. Liischer and P. Weisz, Two loop relation between the bare lattice coupling and the ms 
coupling in pure SU(N) gauge theories, Phys. Lett. B349 (1995) 165-169 
|hep-lat/950200l| . 

[21] M. Liischer, P. Weisz and U. Wolff, A Numerical method to compute the running coupling 
in asymptotically free theories, Nucl. Phys. B359 (1991) 221-243. 

[22] S. Takeda and U. Wolff, Automatic generation of vertices for the Schroedinger functional, 
PoS" LATTICE2007 (2006) 257 [0709 .41671 . 



34 



